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In the heterostructure composed of a topological insulator sandwiched by two s-wave supercon- 
ductors, the time reversal invariant topological superconducting phase, possessing helical Majorana 
edge modes, is found to exist when the two s-wave superconductors form a Josephson junction with a 
7r phase shift. Based on such a heterostructure, a helical Dirac- Majorana interferometer is proposed 
to directly measure a unique transport signature of the helical Majorana modes. Furthermore, we 
envision how our proposal can be realized on the basis of existing materials such as Bi2Se3 or Bi2Te3 
thin films. 
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Recently, due to the discovery of topological insula- 
tors (TIs), the search for new topological phases has at- 
tracted a lot of interests in the condensed matter physics 
community - — . Interestingly, topologically non-trivial 
phases not only exist in insulators, but also in super- 
conductors (SCs), which are dubbed topological super- 
conductors (TSCs). For TIs, the unique feature is the 
existence of topologically protected edge states (or sur- 
face states), and similarly, there are also gapless modes 
at the edge (surface) of the TSCs which are Majorana 
fermions due to particle- hole symmetry. Bound Majo- 
rana modes in a vortex core and chiral Majorana modes 
along the edge were first proposed in p x +ip y SCs 5 - and 
later on helical Majorana (HM) modes have been sug- 
gested to exist in TSCs with time reversal symmetry^"—. 

Majorana fermions may serve as a building block 
for topological quantum computation due to their non- 
Abelian statistics*^. However, unlike TIs, up to now, 
the existence of TSCs and Majorana fermions has not 
been experimentally confirmed. The difficulty of the de- 
tection of Majorana fermions lies in the fact that they 
are charge neutral and have no direct response to electro- 
magnetic fields. Therefore, recently it has been suggested 
to detect Majorana fermions through a Dirac-Majorana 
converter, converting a pair of Majorana fermions into 
one Dirac fermion and vice versaii^. Such a Dirac- 
Majorana converter has been first proposed for the de- 
tection of chiral Majorana modes. However, for the de- 
tection of HM modes, realistic proposals barely exist^. 

In this paper, we consider the heterostructure of a TI 
thin film sandwiched by two SCs. Due to the proximity 
effect on both the top and bottom surfaces of the TI, it is 
found that, if the s-wave pairing functions at the top and 
bottom surfaces have a relative tt phase shift, a topolog- 
ical superconducting phase emerges in the TI thin film 
possessing a HM mode along its edges. Based on such 
a configuration, we propose an experimental setup for a 
HM interferometer and discuss the transport properties 
which can be used to distinguish it from other types of 
Majorana interferometers. 



How can the envisioned heterostructure schematically 
shown in Fig.[TJa) be constructed based on existing mate- 
rials? Superconductivity can be realized by copper doped 
Bi 2 Se<A TlBiTe^, and Bi 2 Te 3 under pressure^, which 
can be straightforwardly integrated into TI thin films 
based on Bi2Se3 or Bi2Te3. Let us describe the mini- 
mal model that captures the essential physics of such a 
system. We assume that the chemical potential lies in 
the bulk gap of the TI, then the low energy physics of 
the TI film is described by two Dirac cones at the top and 
bottom surfaces, given by the four band Hamiltonia n 17 ! 18 

k 

h(k) = mr x + hv f (k x a y r z - k y a x r z ) (1) 

with the field operator ipk = [ci^, c\±, C2^, C2\\ T , where 
1(2) denotes the top (bottom) surface and t (I) denotes 
spin up (spin down); Vf is the Fermi velocity and m 
the hybridization between the top and bottom surface 
states. Here, the Pauli matrices &i denote the spin and 
Ti the opposite surfaces. The Hamiltonian (pQ) satisfies 
time reversal symmetry, THqT~ x = Hq, with the time 
reversal operator T = ia y K (K denotes complex conju- 
gation). The proximity effect from the s-wave SC for the 
top and bottom surfaces can be taken into account by 
the Bogoliubov-de Gennes (BdG) Hamiltonian 



H BdG = Yl *lH BdG (k)* k , * k =( 
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with the s-wave pairing function A given by 
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where Ai( 2 ) is for the top (bottom) surface. The BdG 
Hamiltonian is invariant under charge conjugation C = 



2 



\ X K, where Pauli matrices denote particle and hole. 
To preserve time reversal symmetry, the pairing function 
A 1(2) can only be real, indicating that the relative phase 
of the pairing functions between the top and bottom sur- 
faces can only be or tt. The band dispersion of the 
Hamiltonian (j2j) can be directly solved and we first look 
at some simplified limits. If the chemical potential \i = 0, 
the band dispersion is given by 

E rts = s^Jh 2 v 2 k 2 + E 2 , E g = ^m 2 + A+ + t\ A_ | (4) 

with r, t, s = ± and A± = Al ^ A2 . The above band 
dispersion shows double degeneracy and the two degen- 
erate branches (r = ±) can be traced back to time rever- 
sal symmetry (Kramers partners). We note that for the 
branches with t = +, there is always a gap for any non- 
zero m, A±, while for the branches with t = — , the gap 
will close when m 2 + = A?_. At the critical point, the 
band dispersion ((4]) becomes linear and there are totally 
two Dirac cones (r = ±) which are Kramers partners. 
For one Kramers partner, a single Dirac cone will induce 
the change of the Berry phase by tt between the E g > 
and Eg < regimes, while for the other one, the change 
of the Berry phase is — tt. Therefore, similarly to the 
quantum spin Hall case 19 , a topological phase transition 
is expected to happen across the point E g = in param- 
eter space. When the SC gap A± is much smaller than 
the hybridization gap m, the system should be a trivial 
SC. Thus, we expect that the system should be a topo- 
logical^ non-trivial SC when m 2 + < /S?_ as shown 
in Fig. DJb) (which will be confirmed below). We note 
that the above condition implies if A1A2 < 0. Hence, the 
non-trivial phase only exists when the pairing functions 
of the s-wave SCs for the top and bottom surfaces have 
opposite signs (i.e. they form a Josephson junction with 
a tt phase shift). 

To determine the type of TSC, we consider the limit 
A + = 0, in which we can transform the Hamiltonian 
(0) into the basis c ± , m) = [ci^i) ± c 2,m))> and find 
the obtained Hamiltonian is block diagonal with the two 
blocks related by time reversal symmetry. Each block of 
the new Hamiltonian is exactly the Hamiltonian for the 
chiral TSC discussed, for instance, in Ref. [20|. There- 
fore, our system is expected to be nothing but a TSC 
with HM edge modes. This can be seen in analogy to the 
quantum spin Hall state with helical edge states consist- 
ing of two copies of the chiral quantum anomalous Hall 
stated. To substantiate the above picture, we take into 
account another limit (m = 0) in which the two surfaces 
are decoupled from each other. Then, for one surface, our 
Hamiltonian is the same as the one discussed by Fu and 
Kane in Ref. 21. According to Fu and Kane, a Joseph- 
son junction with a tt phase shift will induce a helical 
Majorana wire along the one-dimensional junction. In 
this decoupled regime (m = 0), the prediction by Fu and 
Kane is consistent with ours that a HM wire emerges at 
the edge of our sandwich structure if there is a tt phase 
shift across the SC-TI-SC Josephson junction. 



Next, we directly solve the Hamiltonian (j2j) in a half 
infinite plane (y < 0) to obtain the effective Hamiltonian 
for the HM mode. We assume that the x-direction is 
translation invariant so that k x is a good quantum num- 
ber. We will first solve the eigen value problem at k x = 

H BdG (k x = 0, -idyMy) = E$(y) (5) 

and then project the Hamiltonian (j2j) onto the subspace 
of the eigenstates at k x = 0. In order to impose open 
boundary conditions, we need to add a quadratic term 
into the hybridization term m = mo + Ed 2 , where B is 
considered to be a small positive number. For simplicity, 
we only consider the case \i = 0, A + = and A_ > 
mo > 0. With such a simplication, we search for the 
zero-energy state (E = 0) with the ansatz <1> oc e Xy (j) for 
the boundary conditions $(0) = and &(y — >• —00) = 0, 
and find the doubly degenerate eigenstates 

*i = 4r ( e ~ Xiy ~ e ~ X2V ) i = ±, (6) 

iVo 

where No is a normalization factor and Ai 5 2 given by 

Ai - 2 = h (~ hVf ± \/ ft2u /- 4B ( mo + A -)) w 

which are real numbers, smaller than 
zero. The wave functions </> + = 

\ [1 + z, 1 - z, -1 - z, 1 - z, 1 - z, 1 + z, -1 + z, 1 + i] T 
and (p- = j [—1 — z, 1 — z, — 1 — z, — 1 + z, — 1 + z, 1 + z, 
— 1 + z, — 1 — z] T are given in the basis e),|l,| 

,e),|2,t,e),|24,e),|l,t,h),|14,h),|2,t,h),|2,;,ft). 
Since the <p± are invariant under charge conjugation 
(C(j)± = (j)±) and can be related to each other by 
time reversal = — </>_ and T</>_ = 0+), they are 

expected to be HM modes. Indeed, for nonzero k x , we 
project the Hamiltonian (j2j) into the subspace of <I>±, 
which can be done by expanding the field operator ^ 
as if(k x ,y) = ^2i®i(y)ji(k x ) (z = ±), and obtain the 
effective Hamiltonian 

^eff = ^HVfk x (^\_(k x )l+(k x ) -7-(M7-(M), ( 8 ) 

k x 

which is exactly the HM mode as expected. 

After determining the condition of the topological 
phase for \i = 0, we can easily extend to the /i ^ regime. 
This is done by analyzing the \i 7^ case and afterwards 
adiabatically connecting the result to the \i = regime. 
By solving the bulk dispersion, the helical topological 
superconductor phase is found to exist in the regime 

\m\ < + (A 2 _-A^). 

In order to confirm the HM mode experimentally, we 
consider the setup shown in Fig. Efa). It is well-known 
that for a TI film, it is possible to obtain the two- 
dimensional quantum spin Hall state with helical edge 
states by tuning the thickness of the film due to quantum 
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FIG. 1: (Color online) (a) Side view of the SC-TI-SC 
heterostructure; (b) phase diagram for the SC-TI-SC het- 
erostructure when the chemical potential /x = 0. 



confinemen t 1 7 i 18 . Then, for the TI film with proper thick- 
ness, we introduce an s-wave superconductor on both the 
top and bottom surfaces near the upper edge as shown by 
the pink regime in Fig. [2fa) (marked by the label TSC). 
By tuning the relative phase of the pairing functions for 
the SCs on the top and bottom surfaces, we can obtain 
a TSC with HM modes at the edges. Hence, at the top 
edge in Fig. Efa), we find one helical Dirac fermion from 
the left hand side (LHS) will be splitted into two Majo- 
rana fermions and then recombined into another helical 
Dirac fermion at the right hand side (RHS). Therefore, 
a HM interferometer is obtained. 

Next, we will discuss the transport properties of the 
HM interferometer within the scattering matrix formal- 
ism. The HM interferometer can be divided into three 
parts, a Y junction at the LHS is connected to another 
Y junction at the RHS by two freely propagating Majo- 
rana wires of different length L\ and L2. The Y junc- 
tion regime consists of one helical Dirac fermion mode 
and two Majorana fermion modes, as shown in Fig.[2fb). 
We name the Dirac fermion mode iV^o)) {ol = e, h) and 
the two Majorana modes |7^ G )) {ol = 1,2), where i (o) 
stands for incoming (outgoing) mode. Time reversal re- 
quires T\i/>f (o) ) = (-)|V£ (i) > and r l7^ o) } = Rl7 a W >> 
while particle- hole symmetry gives C\j™) = I7") and 

C\ifj^ h ^) — IV^^) {v ~ The scattering wave func- 

tion takes the form \^) = ^ ar) {c a ^\^) + d^^)), 
and consequently the S-matrix for the Y junction Sy is 
defined as 



C 2,j,o 
de,j,o 
V dfr j J 



S Y 



( c i,i,i 

d e ,j,i 



(9) 



where j = L,R for the left and right Y junction. With 
the constraints from time reversal symmetry^ 2 -, particle- 
hole symmetry, and the unitary condition, Sy can be 
parametrized as 
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t\. In the 



1*1 1 + 1*2 1 2 + M = 1 5 r l is real and *2 

latter equations, n refers to the reflection amplitude, and 
t\{2) are the transmission amplitudes depending on the 
microscopic details. In the middle regime, the two HM 
fermions are propagating freely, only picking up a phase 
shift, therefore the corresponding S-matrix S m reads 



C2L,i 
ClR,i 
V C 2 RA / 



Sr) 



with 



Tlr 
Trl 



C2L,o 
ClR,o 

\ C 2 Rn / 



j Tlr 



(12) 









J0 2 



(13) 



and Trl = ~T^ R due to time reversal symmetry. The 
angles 61 and 62 denote the phase shifts of the Majorana 
fermions propagating along the two arms of the TSC. 
The phases (in the absence of trapped vortices in the 
TSC region) are given by 0\ = k x L\ + 7r, O2 = k x L2, 
where is the length of the upper (index 1) and the 
lower (index 2) arm in Fig. [2fa). The 7r phase shift, 
which we choose to add to #1, comes from a Berry phase 
of a spin- 1/2 rotation. Combining the above S-matrix for 
three different regions, we can obtain the total S-matrix 
S t defined as 



(14) 



which connects the helical Dirac fermions on the LHS to 
those on the RHS. The matrix elements of St are named 





a(3 



where a, f3 = e,h denote electron or hole and j, / 



L, R denote left or right lead. Then, the conductance is 
given by22r— 



Gil = 



d(Ij 



(15) 



where T' 



■ck/3 



a/3 1 2 



L ji ~ 1 ' re P resen ^ s t ne differential (i.e. 
non-linear) conductance for the current in lead j with 
respect to the voltage in lead I. We only consider the 
zero temperature case for simplicity. Here, TJJ and T^ e 



{j ^ I) are the electron reflection and transmission prob- 
abilities, respectively, while and T?/ 1 {j ^ I) are the 
Andreev reflection and transmission probabilities, respec- 
tively. By solving St from the expressions ([9])- ([T3]h we 



can obtain the T 



af3 



We find T[ e L and T RR vanish be- 



cause of time reversal symmetry. The Andreev reflection, 



4 



Vacuum 




Gll 



FIG. 2: (Color online) (a) Proposed experimental setup 
for the HM interferometer; (b) Zoom into the Y junction, 
where Dirac fermions (full lines) are converted into Major ana 
fermions (dashed lines) and vice versa; (c) Dependence of the 
local conductance Gll — cHl/cIVl and (d) the non-local con- 
ductance Glr — cLIl/cIVr on the phase factors Q 1 / 2 (in units 

of 4). 



the electron transmission, and the Andreev transmission 
probabilities are given by 



rphe rphe 

1 LL — -LRR 



4 r 2 sin 2(*i±f2) 



r\ + 1 - 2r\ cos(0i + 9 2 ) ' 
4M 4 sin 2 (^) 



Tiee rpee "I" 1 ! V 2 

LR — L RL 



r\ + 1 - 2r\ cos(0i + 9 2 ) ' 



1 LR — 1 RL ~ 



r\ + 1 - 2rf cos(0i + 6» 2 ) ' 



(16) 
(17) 
(18) 



respectively, which satisfy the current conservation con- 
dition Y.ip'Ttf = 1- 

From the expressions (fT6]) - ([T8]h we find that the cur- 



rent through our HM interferometer is qualitatively dif- 
ferent from both (i) the current through a TI-SC-TI 
junction 2 - and (ii) the current through a chiral Majorana 
interferomete r 11 ! 12 . First, although there is no electron 
reflection in the HM interferometer, Andreev reflection 
still exists. This is similar to the TI-SC-TI junction^ 
but different to the chiral Majorana interferomete r 11 ! 12 . 
Second, in our case, an electron can transmit both as an 
electron or a hole, which can be used to distinguish the 
HM interferometer from the TI-SC-TI junction 26 , where 
only electron transmission is allowed. 

In the HM interferometer, the normal and Andreev 
transmission probability depend on the phase difference 
between the upper and lower arms. It can be tuned by 
changing the applied bias, by varying the lengths of the 
arms L1/2, or by varying the number of vortices trapped 
in the TSC regime as recently suggested in Ref. 27. The 
dependences of the local conductance Gll and the non- 
local conductance Glr on the parameters 9\ and O2 are 
shown in Fig.[2](c) and (d), respectively. The observation 
of such an interference pattern would be a clear signature 
of HM modes in this setup. 

To summarize, we have proposed a realistic setup for a 
HM interferometer which can be used to confirm the ex- 
istence of HM modes and helical TSCs. The core parts of 
our setup could, for instance, be realized in Bi2Se3 thin 
films selectively doped by copper to turn layers within 
the structure superconducting. A big advantage of our 
proposal is that the superconducting material and the TI 
material are of the same type. Therefore, there should be 
no Schottky barriers between them maximizing the prox- 
imity effect. We have proposed to identify HM modes by 
measuring phase-dependent transport through a HM in- 
terferometer. 
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